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dimensional nonlinear 
Volterra-Fredholm integral equations 
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Abstract 


In this paper, an iterative scheme for extracting approximate solutions 
of two dimensional Volterra-Fredholm integral equations is proposed. Con- 
sidering some conditions on the kernel of the integral equation obtained by 
discretization of the integral equation, the convergence of the approximate 
solution to the exact solution is investigated. Several examples are provided 
to demonstrate the efficiency of the approach. 


Keywords: Volterra-Fredholm integral equation; Iterative method; Dis- 
cretization; Approximation. 


1 Introduction 


The integral equations provide important tools for modeling a wide range 
of phenomena and processes [14], and solving many problems in engineering 
and mechanics which are dependent on finding the solution of their integral 
equations. They are widely used in plasma physics [10], deblurring of two 
dimensional images [8, 20], solving applied boundary value problems [1] and 
Laplace’s equations with boundary conditions [16]. Upon the importance 
of the integral equations, different numerical methods have been developed 
over the years to tackle them, such as time collocation and time discretization 
methods [6, 15], trapezoidal Nystrom method [11], Adomian decomposition 
method [9, 17] and successive iterative scheme [5] but few of them can be used 
for solving two dimensional integral equations such as two dimensional block 
pulse functions [3], finite difference inequalities [19], time-stepping methods 
[7] and block-by-block method [4]. 
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Studies on iterative approaches play an important role to accelerate con- 
vergence rate in solving any system of equations generated by discretizing 
mathematical models in science and engineering problems [5]. The objective 
of this study is to present an iterative approach for extracting approximate 
solutions of two dimensional Volterra-Fredholm integral equations as 


t pb 
u(a,t) = f(a, t) +f :) k(a,t,y, 2, uly, z))dydz, (x,t) € D := [a,b] x [c,d], 
(1) 


where f(z,t) (source function) and k(x,t,y,z,u) (kernel function) are the 
given analytical functions defined on D and D x D x R, respectively. The 
existence and uniqueness of the solution for equation (1) are discussed in 
[12, 15]. This work can be considered as an extension of the method proposed 
in [5]. Note that, the present approach is applicable to a wide class of integral 
equations. The structure of the report is as follows. In Section 2 we transform 
the integral equation into a discretized form. Then, in Section 3, we introduce 
an successful numerical approach which is used subsequently for making up 
the solution algorithm in Section 4. Section 5 demonstrates the efficiency and 
advantages of the proposed algorithm whilst Section 6 concludes the paper. 


2 Integral equation transformation 


Let AQ) _ {a = %0,%1,°°* 5 Un—1, en = b}, A) = {ec = to, t1, aes bis bin = 
d} be equidistance partitions of [a,b] and [c,d], respectively, where hy, = 
Vi41 — Lji,7=0,1,--- n—l and ht =tj41 —t;,7 =0,1,--- ,m—1are the 
discretization parameters of the partitions. Now, if u*(«,t) be an analytical 
solution of (1), then for the partitions AW, A®) on [a,b] and [c, d], we have 


tj b 
wea) = secs i} / Kane ames. O) 


where i = 0,1,--- ,n and j = 0,1,--- ,m. In (2), the integral term can be 
estimated by a numerical method of integration, e.g. Newton-Cotes methods. 
Therefore, by taking equidistance partitions A“), A®), as above with h, = 
Yt — Yi, 1=0,1,--- ,n-lhz = 241-2, § =0,1,--- ,m— 1, and also 
the weights w;, i= 0,1,---,n and wij, r =0,1,--- ,j, equality (2) can be 
written as, 


4, n 
ut; = fig + D> wi wik(2i, tj, yi, zr ut) + O(RY) + O(hH), — (3) 
r=0 I=0 
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where uj, = u*(2i,t;), fig = f(ti,tj), 7=0,1,---,n, 7 =0,1,--- ,m, and 
vy, 4 depends upon the employed method of Newton-Cotes for estimating the 
integral in (2). 

For partitions AW), A®), we consider a nonlinear equations system ob- 
tained by neglecting the truncation error of (2), as follows, 


jon 
, : : 
6,3, =figt w, wik(ez,ts,yp,2%r,€1.r),¢t=0,1,--+,n, 7 =0,1,--+,m, (4) 
oJ oJ Oizo Te 2 ’ 
r=01= 


and suppose that the exact solution of nonlinear system (4) are €7;,7 = 
0,1,--- ,n,j =0,1,---,m. In the following proposition, we seek the condi- 


tions of vanishing |u7; — 7 ;|,¢=0,1,---,n,7 =0,1,--- ,m 


proposition 2.1.Suppose, 
igl> 


(2) | 4 a ql= = man osisn luz 5 — 9 
(ii) k(x, t,y, 2, u(y, z)) € C(D x Dx R), 
k 


(iti) ky (x,t, y,z,u(y, z)) exists on Dx Dx Randy < CHyICaoe where 


= sup |ku(x, t, y, z, u(y, z))|- 


x,ye [a,b] 
t,z€[c,d] 
Then 
ieee |O(hy)| + |O(AZ)| 6) 
pa Spal S755 alld — 0) 


Proof. By (3) and (4), we have 


qn 
= wi(k (Zp, tgs Yls Zr) Usp) — (Xp, tgs Yt» 2rs Ef y)) 
i=0 


r=0 
+O(hy) + O(he). 
According to (iii) 
k(p, ta: YL, Fr, UL y) a k(p, tq; Yl, Zr, fr) 7 Fy (tr tq, YL, Fr, Nr) (Uy ee Gie)s 


(6) 
where for each | = 0,1,--- ,n,r =0,1,--- ,m, m,, is a real number between 
uj, and €/,. Again by (iii) and (6), we conclude that 


qd n 
tna Spal Ss >, S/ wu’, wiluy, — 1+ |O(rg)| + |O(AE)| 
r=0 1=0 
qd n 
< que g — Bgl Y. Dd) wi, wi t+ |O(hY)| + |O(RE)]. 


r=0 1=0 
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Since in every Newton-Cotes formula 77_ Soyo wi, wi = (b— a)(d—c), 


eg lO(ry)| + |ORE)| 


Mea Seal S150) al(d— 0) 


Inequality (5) leads to the following corollary, corollary 2.2. |us ,—& | 
vanishes when h, and h, tend to zero. 


Now, to find the approximate solution, one needs to solve nonlinear equa- 
tion (4). 


3 The successive numerical approach 


Iterative methods are widely used for finding approximate solution of non- 
linear systems of equations [21]. Borzabadi et. al. in [5] presented a succes- 
sive substitution, similar to Gauss-Seidel method for solving one dimensional 
Fredholm integral equation. The nonlinear system of equations (4) has also a 
structure that permits us to approximate its solution by a similar successive 
iterative approach presented in [5]. Hereby we define an iterative process 
which leads to the sequence of matrices {¢(")}. The components of the ma- 
trices satisfy the iteration formula, 


“9: n 
k+1 k 
EO) = fig + > whwik(is ty ws te): (7) 
r=0 1=0 
where 7 = 0,1,--- ,n,j7 = 0,1,---,m and k = 0,1,---. Though, the con- 


vergence scheme can be constructed for detecting approximate solution (4). 
However, we first study the conditions that guarantee the convergence of the 
sequence {€()}, 


theorem 3.1. Considering assumptions of Proposition 2.1, the produced 
sequence {€")} from the iteration process (7) tends to the exact solution of 
(4), say €*, for any arbitrary initial matrix €). 


Proof. By (4) and (7) we have, 
k+1 ace k 
et ) ee = oS So wi wilk(ai, ty, is 26) = k (xi, tj, Yt, Zr, Elin), 


r=0 l=0 


and according to condition (itz) of Proposition 2.1, 
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k+1) 
eee) SS yy: wie = (wa tjs ue 2m (G? - a 
r=0 1=0 
where Th ) is a real number between e) and 7, for 1 = 0,1,--- ,n and 
r =0,1,--:,m. Thus one may obtain the following inequalities 
n 
(k+1 
ef ya igl S O<ien ei? ADS wile 7 (eet ta Yr zr m)| 
0<j<m r=0 1=0 
m n 
xk / 
sy ofien ay =. 5 SoS owhwi, 
0<j<m r=0 1=0 


where 7 = 0,1,--- ,n,j = 0,1,---,m. By setting A = 7(b — a)(d —c) we 
conclude that 


(k+1) * * 
aie Fl. 
an ees. hes ee) - & 
0<j<m 1<j<m 
By mathematical induction on k, we get 
(k+1) _ k * 
max <\* max ie 
0<i<n IS 5 va O<i<n ee ~ Sigh; 
0<j<m 0<j<m 
for each k = 0,1,---. Since 0 < A < 1, then, k 4 +00 implies that 
max o<i<n |&; ; (+1) — €};| vanishes. 
0<j<m ; 


4 Algorithm of the approach 


In this section, we propose an algorithm on the basis of the above discussions 
to solve the Volterra-Fredholm integral equation (1). This algorithm is pre- 
sented in two stages, the initialization and the main steps. 


Initialization 

Choose € > 0, and equidistance partitions AY) = {a = x = yo,21 = 
Yis*+ 5 @n-1 = Yn-1,Ln = Yn = 0} on [a,b] with the step size hy = 
T4+1 — VM, 7=0,1,--- yn—-1, plus A(?) {c to Zo, t1 A1,°** bt: = 
Zn—1,tn = 2n = d} on [c,d] with the step size hy = tj41 —tj), j = 
0,1,---,m—1 and an initial matrix €. Set k = 0 and go to the main 
steps. 


Main steps 
Step 1. Compute é“*+1) by (6), and go to Step 2. 
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Step 2. Compute max o<i<n jae) - ae 


aj 
O0<j<m 


Step 3. If Max 0<i<n ee _ 
O<j<m 


| and go to Step 3. 


(| < ¢€, stop; Otherwise, set k =k +1 and 
go to step l. 


In the next section the advantages and the influence of the proposed 
approach in thrilling convergence rate of the solution for the problems is 
demonstrated via some examples. 


5 Numerical examples 


Suppose u*(a,t) is the exact solution of Volterra-Fredholm integral equation 
(1) and &s5; i=0,1,---,n, 7 =0,1,--- ,misasolution obtained by applying 
the given algorithm with a known ¢ > 0 and partitions A™ and A®). To 
compare the precision of the approximate solution, the discrete error function 


ela, t;) = lear eats) ~ €(ai, tj)], 1=0,1,-:- »n, j =0,1,--- »™, (8) 
is established. 


Example 5.1. In this example, we apply the developed method to a two 
dimensional Fredholm integral equation as follows [13], 


1 x BAR <p 
—(l+y+z)u?(y, z)dydz. 
(Gt+z+t? 6a+%) if | pee tT UE Wey 
This integral equation has analytical solution u(x,t) = eens on [0, 1] x 
[0,1]. We take « = 10~® and partitions with the discretization parameters 
hy = zag and hy = 74g. The initial matrix €) = 0 is considered first to 
start the algorithm. In Table 1, one can see all acceptable values for error 


estimation (8) which is obtained by applying the developed algorithm to the 
illustrated equation. 


u(a,t) = 


Table 1: Error estimation for Example 5.1 
at 0.0 0.2 0.4 0.6 0.8 1.0 


0.0 0.000 0.000 0.000 0.000 0.000 0.000 

0.2 | 3.980x 10-§ 3.317x 107-8 2.843x 10-8 2.488 10-& 2.211107 1.990 x 107° 
0.4 | 7.958 10-& 6.631 x 107-8 5.684x 10-8 4973x1078 4.421 10-§ 30979 x 107-6 
0.6 | 1.193 x 107° 9.943 x 10-& 8.522x10-® 7.457x 107° 6.629x 107° 5.966 x 107 
0.8 | 1.590 10-® 1.325 1075 1.136x 10-5 9.940x 10-& 8.835 x 10-® 7.952 x 10-8 
1.0 | 1.987x 1075 1.656x10-® 1.420x 1075 1.242 1075 1104x1075 9.937 x 1076 


Example 5.2. In this example, we apply our method for the following two 
dimensional Fredholm integral equation [2], 
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u(x,t) = 2e* — a — 75" + zr +f | (xy + te*)u(y, z)dydz. 
0 Jo 


The analytical solution of this integral equation is u(z,t) = ve~* + t 
on [0,1] x [0,1]. By solving this equation, we observe that the proposed 
algorithm does not give rise to a convergent sequence. So, to overcome this 
shortcoming, we put [0,0.1] x [0,0.1] in place of [0,1] x [0,1], where the 
conditions of Theorem 3.1 hold. Then 


, , 1799 43 9 
t 


44 got 4 t -1 
2000 120000° ' 100 


t) = Ze- ee 
u(x,t) = xe 3000 7° 


0.1 p01 
+f | (xy + te*)u(y, z)dydz. 
0 Jo 


Table 2 shows that in this region of integration, approximate solution tracks 
the exact one, almost precise. 


Table 2: Error estimation for Example 5.2 
at 0.0 0.02 0.04 0.06 0.08 0.1 


0.0 0.000 3.598 x 10- 7.197 x 10— 1.080 x 10~ 1.439 x 10- 1.799 x 10- 

0.02 | 3.178 x 10-!! 6.776 x 10-!! 1.037 x 10-19 -1.397 x 107-19 -1..757 x 107! 2.117 x 10710 
0.04 | 6.356 x 107!!! 9.954 x 107" 1.355 107! 1.715 x 107! 2075 x 10-19 =. 2435 x 1071° 
0.06 | 9.533 x 107! 1.313 x 107-19 -1.673 x 10-19 2.033 x 10-19 2.393 x 10-!9 2.752 x 10710 
0.08 | 1.271 x 10-1 1.631 x 10-1 1.991 x 107! 2.351 x 10-!° 2.710 x 10-!° 3.070 x 10~!° 
0.1 | 1.589 x 10-9 1.949 x 10-1 -2.308 x 10-1 2.668 x 10-1 3.028 x 10-19 3.388 x 10-10 


Example 5.3. In this example, we apply the proposed method to the fol- 
lowing two dimensional Volterra integral equation [18], 


z*sin?(t), 2° ,sin(2t) 


u(x,t) = xsin(t)(1 9 )+ 10! 5 t) 


+f fev + goal 2) u2(y, 2)dyde. 


This integral equation has analytical solution u(x,t) = xsin(t) on [0, 1] x 
[0,1]. We take « = 10~® and partitions with the discretization parameters 
hy = za and hy = sav The initial matrix €() = 0 is considered for starting 
the algorithm. Table 3 illustrates the precision of the approximate solution 
by showing the error criteria (7) corresponding to the given partition. 


Table 3: Error estimation for Example 5.3 

at 0.0 0.2 0.4 0.6 0.8 1.0 

0.0 | 3.906 x 10-8! 2.651 x 10-8 —-1.716x 10°" ~=7.961 x 10°? 2.432 x 10-7 4.102 x 10-* 
0.2 | 3.227 x 10788) 1.722 1078 ~— 8.025 x 1078 )=— 2.172 x 1077 4.268 x 1077 «6.769 x 1077 
0.4 | 4.711 x 10788) 8.889 x 1078 ~—- 2.616 x 1077 —5.724 x 1077 1.025 x 1078 ~—-1.569 x 1078 
0.6 | 1.199 x 10782 2.928 x 10-7 = 7.087 x 1077 —-1.838 x 1076 =. 2.213 x 1078 ~—- 3.2938 x 1078 
0.8 | 5.513 x 10-82, 7.761 x 10-7 —-1.742x 10-© = 3.063 x 10-6 3=— 4.852 x 10-6 )~—7.135 x 10-8 
1.0 | 3.906 x 10-3! 1.833 x 10-® ~=— 4.003 x 10-© ~— 6.847 x 10 }~=—- 1.067 x 10-5 1.571 x 107° 
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Example 5.4. In this example, we apply our method to a Volterra-Fredholm 


integral equation as follows [3], 


1 1 t 1 
u(a,t) =a? +<at at? zt? + | aty?z? u(y, z)dydz. 
15 16 o Jo 


This integral equation has analytical solution u(«,t) = 2?+<2t on [0,1] x [0, 1]. 
We take e = 10~® and partitions with the discretization parameters hy = na 
and hy = jaa" The initial matrix € = 0 is considered for starting the 
algorithm. Table 4 exhibits good error values by applying the developed 
algorithm. 


Table 4: Error estimation for Example 5.4 
xt! 0.0 0.2 0.4 0.6 0.8 1.0 


0.0 | 0.000 0.000 0.000 0.000 0.000 0.000 

0.2} 0.000 4.063 10-8 2.567x 10-7 9.000x10-7 2.425x10-® 5.598 x 10-6 
0.4} 0.000 8126 10-8 5.133x 1077 1.800x10-® 4.850x10-® 1.119 x 10-> 
0.6 | 0.000 1.219 10-7 7.699x 107-7 2.700x10-® 7.276107" 1.679 x 107° 
0.8 | 0.000 1.625 10-7 1.027x10-® 3.600x10-® 9.701 x 107° 2.239 x 107° 
1.0 | 0.000 2.031 x 107-7 1.283x107-® 4.500x10-® 1.213x 1075 2.799 x 10-® 


6 Conclusions 


In this paper, an iterative approach for obtaining approximate solutions for 
two dimensional Volterra-Fredholm integral equations, considering some spe- 
cial conditions on the kernel, as continuous differentiability of kernel, is pro- 
posed. Theorem 3.1 provides a sufficient condition for convergence of the 
approach, but it is not necessary. Therefore, Examples 5.1, 5.3, and 5.4 show 
that, despite the lack of conditions, convergence of the proposed method 
holds for a class of two dimensional Volterra-Fredholm integral equations. 
Also the changing in problem for holding conditions of Theorem 3.1 lead to 
the convergence of the method, as it is described in Example 5.2. The valid- 
ity and efficiency of the proposed scheme is demonstrated on the examples 
included. 
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